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The library for Mathworks® Simscape™ does not currently contain a model for a closed- 
volume fluid tank where the ullage pressure is variable. In order to model a closed-volume 
variable ullage pressure tank, it was necessary to consider at least two separate cases: a 
vertical cylinder, and a sphere. Using library components, it was possible to construct a rough 
model for the cylindrical tank. It was not possible to construct a model for a spherical tank, 
using library components, due to the variable area. It was decided that, for these cases, it 
would be preferable to create a custom library component to represent each case, using the 
Simscape language. Once completed, the components were added to models, where filling and 
draining the tanks could be simulated. When the models were performing as expected, it was 
necessary to generate code from the models and run them in Trick (a real-time simulation 
program). The data output from Trick was then compared to the output from Simscape and 
found to be within acceptable limits. 
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I. Introduction 

M athWorks® Simscape™ (Simscape), is modeling and simulation software that can be used to represent a variety 
of physical systems. The National Aeronautic and Space Administration (NASA) is investigating further use of 
Simscape by the Application & Simulation Software Engineering Branch, NE-C1, at Kennedy Space Center (KSC). 
Simscape models must be proven to be capable of running in real-time, using Trick (a real-time simulation program), 
and with an acceptable level of accuracy. I have been tasked with creating a functional model of a fluid tank, where 
the gas pressure inside of the tank will increase or decrease with a corresponding change in liquid volume. There are 
two major tank geometries that will be considered for this model: a vertical cylinder, and a sphere. 


II. Simscape Components and Models 


A. Requirements 

1. The component should have one liquid inlet/outlet, located at the bottom. 

2. The tank pressure should vary with changes in the liquid volume. 

3. The component should account for heat transfer between the tank and the liquid in the tank. 

4. The component will be considered to be well insulated, so the heat transfer effects with the environment will 
be assumed to be negligible. 


B. Rough Model 

In order to have a basis for 
comparison, I constructed a model 
that is representative of a 
cylindrical tank, using Simscape 
library components. It uses a 
hydraulic cylinder that is physically 
connected to a pneumatic cylinder 
through the piston rods of the 
cylinders. When using this model, if 
the dimensions of the tank need to 
be changed, the changes have to be 
incorporated into both cylinders, 
this is more time consuming and 
could introduce the possibility for 
errors. Additionally, the thermal 
port on the pneumatic cylinder is 
meant to interact with the gas, and I 
am more interested in the heat 
transfer between the liquid and the 
tank. 



Figure 1. Image of rough Simscape model. 
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C. Custom Cylindrical Tank Component 

Using the Simscape language, it is possible to modify some library components, or, to 
create complete custom components. In order to satisfy the requirements for the fluid tank, 
I decided to start from the beginning. To simplify the equations, I made certain reasonable 
assumptions: The gas in the tank is assumed to be an ideal gas, temperature change of the 
gas due to volume change (compression or expansion), is considered to be negligible, heat 
transfer between the gas and liquid is considered to be negligible, the tank is considered to 
be well insulated, change of tank volume due to pressure change is not considered. The 
port labels for the component require some clarification. H is the hydraulic port, where 
the fluid flows into or out of the tank. A is a physical signal port that outputs the contact 
surface area between the tank and the fluid, it is used with the custom heat transfer 
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Figure 2. Image of custom 
Simscape cylindrical tank 
component. 
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components during simulation. V is a physical signal port that outputs the fluid volume. T is a thermal port that is 
used to model the heat transfer, between the tank and the fluid in the tank. P is a physical signal port that outputs the 
pressure in the tank. M is a physical signal port that outputs the fluid mass, and it is used with the custom thermal 
mass component to model heat transfer. Temp is a physical signal port that outputs the fluid temperature. 


D. Cylindrical Tank Equations 

The following equations were used to create the cylindrical tank component. 
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E. Testing Custom Tank Models 

In order to verify the functionality of 
the custom tank components, they were 
added to a model, where the interaction 
with other components could be tested 
and verified. In the model, the “Variable 
Head Tank” supplies the fluid, and the 
“Hydraulic Constant Pressure Source” 
provides the pressure to initiate flow. 
Heat transfer in the tank component is 
modeled through the thermal port “T”, the 
thermal mass associated with the fluid in 
the tank is represented by “Thermal 
Mass”, and the thermal mass of the tank is 
represented by “Tank Mass”. The heat 
transfer between the tank and the fluid is 
represented by the “Convective Heat 
Transfer” and “Conductive Heat 
Transfer” components. 

□ AAA 


F. Custom Spherical Tank Component 

The assumptions for the spherical tank are the same as the assumptions listed for the 
cylindrical tank. In order to represent the relationship between the area change and the 
height change of the fluid, it was necessary to develop this component from the 
beginning, using the Simscape language. The inputs and outputs for this tank are similar 
to those found on the cylindrical tank, described in section (II. C). 



Figure 4. Image of custom spherical 
tank Simscape component. 
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G. Spherical Tank Equations 

The following equations were used to create the spherical tank component. The derivations for the volume and 
area equations are in the Appendix. 


A = n(2rh f - h f 2 ) 
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H. Custom Library 

After writing the code for 
components in Simscape language, it is 
necessary to save them in a separate folder. 

In order to be considered as a Simscape 
library, this folder must begin with a “+”, so, 
an example would be: +Folder_name, I 
would suggest adding this directory to the 
Matlab path. From the Matlab command 
line, the command to compile and build the 
components in the library folder is: ssc_build 
Folder_name. If there are no errors in the 
Simscape code files, the library should 
resemble Fig. 3. The default image for a 
custom Simscape component is a labeled box 
with the corresponding inputs and outputs. If 
a different image is desirable, it is relatively 
simple to add. Save a copy of the image in 
the “+Folder_name” directory. The file name 
for the image has to match the file name of 
the component, for example: for the 

Simscape file, “custom file.ssc” the image would need to be “custom file.img”, where img is a supported file format 
for the image. The supported image file formats are: jpg, bmp, and png 3 . 
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Figure 5. Image of custom component library. 
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I. Additional Custom Components 

When attempting to model the heat transfer between the fluid and the tank walls, I found that the standard library 
components were missing some necessary functionality. In the “Thermal Mass” standard library component, the mass 
is a fixed parameter that cannot easily be adjusted during simulation. For the tank components, the fluid mass is 
changing continuously with fluid volume. In order to consider the change in mass, I modified the standard “Thermal 
Mass” component to accept the value for mass as a physical signal input. Also, both tank components were modified 
to output the fluid mass as a physical signal. A similar problem was encountered with the heat transfer library 
components (Conductive Heat Transfer and Convective Heat Transfer), where the area is a fixed parameter. The 
surface area of the liquid that contacts the tank changes continuously with fluid volume as well. A similar solution 
was employed to account for the changing area, the heat transfer components were modified to accept the area as a 
physical signal input, and the tank components were modified to output the area as a physical signal. 

III. Verification of Models 


A. Cylindrical Tank Comparison 

The cylindrical tank component was incorporated into a model, similar to the model shown in Fig. 3. The 
simulation was run, and the results were recorded. This procedure was repeated for the rough model, shown in Fig. 1. 


Next, some of the results from the two simulations were plotted and compared, shown in Fig. 6 and 7. 



I then compared the two simulations, using the “Simulation Data Inspector” in Simscape, the results are shown in Fig. 
8. The two runs match each other closely, the Tank Volume curves are within 0.09%, the Pressure curves are within 
2.26%, and the Flow curves are within 1.78%. 


B. Code Generation 

After comparing the two models and finding the results to be satisfactory, the next step is to generate code from 
the Simscape model, so it can be run in the Trick real-time simulation environment. The results from the Trick 
simulation will then be compared to the results from the Simscape simulation, to determine if any variation exists 
between the two simulations. Also, it is important that the model is “real-time capable”, meaning that the computations 
performed during one second of simulation time must take less than one second of real-time to complete. 

Certain configuration parameters in Simscape must be set before code generation can be done. I saved a copy of 
the model onto the “KSC Template”, where the parameters were already set to the required values. I then ran the 
simulation to verify that it would 
compile and run with the new settings. 

In the Simscape model, I replaced all 
of the Scopes with output ports, and 
used Pulse Generators to generate an 
input signal. Using Matlab Coder, as a 
test, I was able to successfully generate 
C++ code from the model. In order to 
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generate all of the necessary files to run the simulation in Trick, a series of steps must be taken. First, I ran 
sim_mc_buildScript with the appropriate command line arguments and file paths, from a terminal. If there were any 
errors in the build, I addressed them, and repeated the command until the build was successful. Next, I created a 
data_logging.txt file with input and output names that matched the model inputs and outputs, using the 
SubsystemTestBench Excel template. I then saved a copy of the data_logging.txt file into file_path/SIM/data. Next, 
I included the data_logging.txt file in the header of file jath/SIM/RUN/input. I ran the simtcidbuildScript with the 
appropriate command line arguments from a terminal, if any errors appeared, I addressed the errors and repeated until 
the script ran successfully. I verified that an S_main_operatingsystem.exe file existed in the file_path/SIM directory. 

To run the simulation in Trick, I changed directory to file_path/SIM. The next step was to run the 
S_main_operatingsystem.exe file 
with argument RUN/input. Using 
Ctrl Z to stop the command, and 
then typed “bg” (to move the 
program to the background). I 
opened the Trick simulation control 
by typing: sim_control.tcl -m 

machine host name -port number. 

The Trick control panel opened 
with the simulation paused. I ran the 
simulation for a sufficient time to 
gather the necessary data. Once the 
simulation was complete, there was 
a log_filename.trk file in the 
file_path/SIM/RUN directory. I 
converted this file to a comma- 
separated value (CSV) file, using 
the command trk2csv with the 
appropriate arguments 4 . Next, I 
imported the CSV data and 
compared it to data that was saved 
from a Simscape simulation, using Fi § ure 9 - Comparison of Simscape data and Trick data. 

Matlab to plot the results. 

As shown in Fig. 9, the data from the Trick simulation overlays the plot from the Simscape simulation exactly. 
The results shown are from the cylindrical tank model, the results from the spherical tank model were similar when 
the simulations were compared. 


Flow Rate 



IV. Conclusion 

When comparing the custom cylindrical tank to the one provided Simscape library components, the results are 
promising. Unfortunately, I have no rough Simscape model to compare the custom spherical tank model against, to 
verify the accuracy of the model, other than the equations that define the tank geometry, though the model performs 
as expected. When the need arises, The use of the Simscape language to create custom components, or to extend 
existing Simscape library components is straight forward, feasible and recommended. 
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Appendix 

Derivation of Spherical Tank equations: 

r(/i) 1 2 * 4 = R 2 - (R - /i) 2 , 

Or 

r(/i) 2 = 2 Rh — h 2 

Where r(h) is the radius of the fluid level at height h and R is the tank radius. 
So, the volume would be: 

V = J A dh, 

Where A = n * r(/i) 2 

rh r Hq 

V = I (2nRh - nh 2 ) dh + I (2nRh - nh 2 ) dh 

Jo 

V = n^R{h 2 - hi) - - — + Rhl - y| 
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